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Abstract 

We  report  two  photon  laser  induced  fluorescence  measurements  near  the  nozzle  exit  and  in  the  plume  of  a 
1  kW  arcjet  thruster.  Atomic  hydrogen  number  densities,  axial  and  radial  velocities,  and  translational 
temperatures  are  measured  in  the  expansion  plume  of  a  1  kW  arcjet  operating  on- hydrogen  propellant. 
Data  from  recent  computational  models  are  compared  with  the  atomic  hydrogen  density  measurements. 


Introduction 

Arcjets  are  expected  to  play  an  ever  increasing  role 
in  satellite  propulsion  needs,  primarily  for 
stationkeeping  and  on-orbit  maneuvering  in  the  near 
term.  While  the  technology  is  considered  viable 
enough  to  be  deployed  on  a  Telstar  IV 
communications  satellite  for  stationkeeping,1  arcjet 
technology  is  far  from  maturity.  In  order  to  compete 
successfully  with  chemical  propulsion  systems  for 
on-orbit  missions,  further  improvements  in  arcjet 
propulsion  systems  are  still  required."  If  needed 
improvements  in  the  performance  level  and 
efficiency  of  arcjets  are  to  be  achieved,  an  increased 
understanding  of  the  fundamental  physical  processes 
that  govern  the  operation  of  an  arcjet  is  essential.  In 
addition,  an  ability  to  predict  the  plume  behavior  of 
a  space  propulsion  device  is  necessary  for  prediction 
and  amelioration  of  spacecraft  interactions. 
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Significant  arcjet  energy  loss  results  from  velocity 
profile  losses  due  to  thick  internal  boundary  layers 
in  the  arcjet  nozzle  and  from  frozen  flow  losses  such 
as  molecular  dissociation.  To  quantify  profile 
losses,  both  gas  velocity  and  density  distributions 
must  be  known.  In  addition,  arcjet  models,  which 
are  necessary  for  timely  and  cost-effective 
improvements  to  design,  must  be  tested  by 
comparison  with  key  physical  parameters. 

At  present,  only  limited  inroads  have  been  made 
into  the  problem  of  plume  density  measurements. 
For  hydrogen  arcjet  thrusters,  determination  of 
species  density  has  previously  only  been 
accomplished  through  a  few  diagnostic  techniques. 

Our  previous  work  using  multi-photon  laser  induced 
fluorescence  (LIF)  provided  the  first  measurement 
of  ground  state  hydrogen  in  an  arcjet  plume.3 
Atomic  hydrogen  density  profiles  were  reported 
near  nozzle  exit  and  a  preliminary  calibration  was 
performed  to  convert  the  relative  density  profiles  to 
absolute  number  densities. 

Knowledge  of  the  atomic  density  profile  at  nozzle 
exit  is  important  for  evaluating  our  understanding  of 
the  physical  processes  in  the  arciet  thruster  and  the 
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computational  models  that  have  been  written  to 
describe  these  processes.4*9  The  density  of  atoms 
relative  to  the  density  of  molecules  is  an  indication 
of  how  much  energy  is  lost  into  dissociation  of  the 
hydrogen  molecules  and  not  recovered  through 
recombination  into  translational  kinetic  modes. 

To  determine  the  molecular  dissociation  fraction, 
the  molecular  species  density  needs  to  be  known  as 
well  as  the  atomic  number  density.  The  recent  use 
of  Raman  spectroscopy10,11  in  an  arcjet  plume  has 
provided  the  first  information  on  molecular 
hydrogen  densities  at  the  arcjet  nozzle  plane. 

In  addition  to  LIF  on  the  ground  state  atoms,  other 
absorption  spectroscopy  techniques12'13  have  been 
used  on  arcjet  plumes  to  characterize  hydrogen  atom 
densities.  These  VUV  and  XUV  spectroscopy 
approaches  are  quite  difficult  to  implement  in 
practice  and  are  limited  to  determination  of  line-of- 
sight  averaged  number  densities  at  downstream 
locations  in  the  plume  where  the  optical  depth  is  not 
too  high. 

For  measuring  velocity  and  temperature,  both, 
ground  state  LEF3  and  excited  state  LEF  4'14'13  have 
proven  to  be  accurate  and  essentially  nonintrusive. 
Though  the  excited  states  of  hydrogen  are  more 
accessible,  most  atoms  in  the  plume  region  are 
expected  to  be  in  the  ground  state.  These  two  types 
of  LIF  temperature  and  velocity  measurements  allow 
examination  of  the  differences  between  the  excited 
state  and  ground  state  species  in  the  non-equilibrium 
plume  environment. 

Experiment 

This  work  is  performed  using  two-photon  LIF 
(2PLIF)  on  ground-state  hydrogen  atoms  in  the 
plume  of  a  1-kW-class  hydrogen  arcjet.  Since  any 
propellant  gas  is  opaque  for  the  VUV  radiation 
required  for  single-photon  excitation  of  ground  state 
atoms,  2PLIF  at  correspondingly  longer  wavelengths 
is  preferred.  The  diagnostic  technique,  which  has 
been  developed  primarily  for  use  in  flames,16'17'18 
uses  two  photons  at  205  nm  to  promote  the  atoms 
from  the  n=l  to  the  n=3  electronic  state  (the  L{3 
transition).  Subsequent  3-2  fluorescence  is  observed 
(the  Ha  transition),  as  indicated  in  the  energy  level 
diagram  in  Figure  1 . 
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The  apparatus  is  shown  schematically  in  Figure  2. 
The  laser  is  a  pulsed  dye  laser  ,  pumped  by  a 
Nd:YAG  with  a  repetition  rate  of  10  Hz  and  a  pulse 
width  of  6  ns.  The  dye  laser  output  at  615  nm  is 
frequency  tripled  to  achieve  about  0.5  mJ  per  pulse 
at  205  nm.  A  mirror  turns  about  80%  of  the  beam 
toward  the  arcjet  chamber  through  a  variable 
attenuator  placed  in  the  beam,  the  remaining  20%  of 
beam  energy  is  directed  toward  a  microwave- 
discharge  source  of  atomic  hydrogen'. 

For  axial  velocity  measurements,  the  beam  is  sent 
directly  down  the  axis  of  the  arcjet  flow  (Path  1)  and 
is  focused  with  a  telescope  lens  configuration  (not 
shown)  outside  the  chamber  with  a  focal  length  on 
the  order  of  2  m.  For  radial  measurements,  the 
unfocussed  beam  is  sent  to  a  turning  prism  inside 
the  chamber  located  underneath  the  arcjet  (Path  2), 
directed  to  pass  vertically  through  the  plume,  and  is 
focused  by  a  200  mm  lens.  The  laser  beam  and 
optics  remain  fixed,  while  the  arcjet  is  mounted  on  a 
motion  control  x,  y,  z  stage  to  translate  it  for  probing 
different  regions  of  the  plume.  . 

A  filtered  photomultiplier  tube  (PMT)  is  placed 
behind  the  final  turning  mirror  before  the  chamber 
in  order  to  detect  amplified  spontaneous  emission 
(ASE)  that  may  propagate  back  along  the  laser  beam 
path  (see  further  discussion  below).  A  200  mm 
focal  length,  2"  diameter  lens  is  placed  inside  the 
chamber  to  collimate  the  LIF  that  is  emitted  toward 
the  side  window.  The  light  is  collected  outside  the 
chamber,  focused  through  a  1  mm  aperture,  and 
detected  with  a  filtered,  gated  PMT.  Since  the  LIF 
occurs  at  656  nm,  the  filters  used  are  a  656  nm 
bandpass  interference  filter  and  an  RG  645  color 
glass  filter  (thus  all  scattered  laser  light  (205nm)  is 
filtered  out).  The  gated  PMT  is  an  ordinary 
Hamamatsu  928  tube  with  a  special  socket  that  is 
triggered  to  detect  light  for  2  (is  during  each  laser 
pulse  and  is  off  between  pulses.  This  is  an 
important  feature  that  allows  us  to  operate  the  PMT 
at  full  voltage  without  exceeding  the  anode  current 
limit  due  to  the  bright  arcjet  emission  background. 
Saturation  is  avoided  by  setting  the  voltage  level  so 
that  both  average  current  and  peak  current  during 
the  gate  remain  below  the  maximum  specified  for 
the  PMT.  A  gated  integrator  with  a  30  ns  gate  is 
used  to  amplify  and  average  the  Ha  fluorescence 
seen  by  the  PMT.  Alternatively,  the  PMT  signal  can 
be  digitized  by  a  fast  oscilloscope  to  obtain 
fluorescence  lifetimes  and  quenching  information. 


The  weak  UV  beam  that  is  sent  into  the  discharge 
cell  is  focused  with  a  150  mm  lens.  The  cell  is  run 
with  a  slow  flow  of  a  few  Torr  of  helium  carrier  gas 
and  a  few  percent  hydrogen.  The  LEF  is  detected 
through  a  filtered  (ungated)  PMT.  Simultaneous 
detection  of  LIF  from  the  cell  during  each  spectral 
scan  of  the  arcjet  LIF  provides  a  zero-velocity 
comparison  from  which  to  measure  Doppler  shifts. 
This  discharge  cell  provides  wavelength  calibration 
but  not  density  calibration.  A  separate  discharge 
cell  is  used  for  density  calibration  inside  the 
chamber  when  the  arcjet  is  not  operating. 

9 

The  procedure  for  obtaining  absolute  number 
density  is  based  on  a  method  reported  by  Meier  et 
al.  '  •  After  a  relative  number  density  scan  and 

corresponding  lifetime  data  are  taken,  the  arcjet  is 
turned  off  and  translated  away  from  the  collection 
volume,  the  chamber  is  opened,  and  a  second 
calibration  cell  is  placed  in  the  detection  volume 
with  the  laser  beam  passing  through  it.  The  LIF 
signal  and  lifetime  is  measured  for  the  hydrogen 
atoms  present  in  the  cell.  Atoms  are  present  due  to 
the  flow  of  hydrogen  gas  through  a  microwave 
discharge  prior  to  entering  the  cell.  Since  the  laser 
beam,  detection  optics  and  electronics  have 
remained  the  same,  the  LIF  from  the  arcjet  and  from 
the  cell  will  have  the  same  proportionality  to 
absolute  number  density  after  correcting  for  any 
differences  in  quenching: 

qA 

Na - •  Nc  -C 

iN  H  —  gC  H 

A  c 

S  and  S  are  the  signal  sizes  (integrated  over  the 

entire  spectral  width)  from  the  arcjet  and  the  cell, 
a  c 

Nh  and  Nh  are  the  absolute  number  densities  of 
atomic  hydrogen  in  the  arcjet  and  the  cell,  and  CQ  is 
the  scaling  factor  for  the  difference  in  quenching. 
CQ  is  given  by  the  ratio  of  fluorescence  decay  rate  in 
the  arcjet  to  that  in  the  calibration  cell. 

The  absolute  number  density  of  hydrogen  atoms  in 
the  discharge  cell  is  obtained  using  a  standard 
chemical  titration  method.  A  schematic  of  the 
calibration  Cell  used  for  this  method  is  shown  in 
Figure  3.  The  cell  is  made  out  of  glass  tubing  and  is 
approximately  30  cm  (12  in.)  long  and  5  cm  (2  in.) 
in  diameter.  Important  to  the  design  of  the  cell  is 
the  long  drift  tube  area  with  a  Teflon  tubing  liner  to 
reduce  gas/wall-  interactions.  This  drift  tube  allows 
complete  mixing  of  the  titration  gas  and  the  gas  of 


interest.  In  this  case,  the  titration  gas  is  NO?  and  H 
is  the  gas  to  be  calibrated. 

The  hydrogen  (with  a  helium  carrier  gas)  enters  the 
drift  tube  area  near  one  end  after  passing  through  a 
microwave  discharge.  A  dilute  mixture  of  NO?  in 
helium  is  added  to  the  flow  through  a  long  thin  tube 
and  enters  the  larger  drift  tube  region  at  a  location 
past  the  point  where  the  hydrogen  has  entered.  This 
reduces  the  amount  of  NO?  that  might  diffuse 
upstream  and  enter  the  microwave  discharge. 

The  N02  reacts  rapidly  with  hydrogen  atoms  when 
the  two  are  completely  mixed  in  the  large  drift  tube 
volume.  The  calibration  involved  flowing  1  slpm  of 
a  2%  hydrogen  in  helium  mixture  and  then  the 
addition  of  2%  N02  in  helium  until  the  signal 
decreased,  typically  at  about  .08  slpm  of  N02/He.  A 
small  vacuum  pump  brought  the  cell  pressure  down 
to  950  Pa  (7.2  Torr)  during  gas  flow. 

The  LEF  signal  decreases  until  it  is  gone  at  the  point 
where  the  partial  pressure  of  added  N02  is  equal  to 
the  partial  pressure  of  hydrogen  atoms  in  the  cell. 
Control  of  the  amount  of  N02  added  allows  the 
determination  of  the  amount  of  hydrogen  atoms 
present  in  the  calibration  cell.  Figure  4  shows  data 
from  one  calibration  run.  The  NO?  concentration  at 
the  x  intercept  is  equal  to  the  concentration  of 
atomic  hydrogen  present  as  each  hydrogen  atom  is 
reacted  away  for  each  NO?  molecule  added  in  the 
fast  reaction: 

H+NO?  -4  OH+NO 

The  large  amount  of  buffer  gas  and  low  amount  of  H 
significantly  reduce  possible  secondary  reactions, 
making  them  negligible  for  the  purposes  of  this 
experiment.  Calibration  took  place  either  before  or 
after  each  arcjet  firing,  usually  on  the  same  day. 

All  thruster  data  were  taken  on  a  1-kW-class  arcjet 
designed  by  NASA  Lewis  Research  Center.21  The 
cathode  gap  was  set  to  0.07”.  The  arcjet  operated  at 
137 V,  10A  (1.37  kW)  on  a  hydrogen  gas  flow  of 
13.1  mg/s  (8.74  slpm)  at  an  operating  chamber 
pressure  of  6  Pa  (45  mTorr).  The  gas  flow  and 
power  were  chosen  to  closely  match  the  conditions 
under  which  the  most  complete  set  of  previous 
diagnostic  data  for  a  1  kW  hydrogen  arcjet  were 
taken,  in  order  to  facilitate  comparisons  between 
different  data  sets  and  between  data  and 
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models.10'15'22  A  voltage-current  curve  of  our  arcjet, 
was  shown  in  a  previous  publication/  The 
hydrogen  flow  controller  for  the  arcjet  as  well  as  the 
flow  meters  for  the  discharge  cell  were  calibrated 
using  a  wet  test  meter. 

Results -and  Discussion 

Figure  5  shows  a  sample  spectral  scan  of  the 
hydrogen  atom  profile  from  the  arcjet  and  from  the 
discharge  cell,  taken  with  an  axial  laser  beam.  Each 
profile  is  fit  to  a  Gaussian  shape  using  a  Levenberg- 
Marquardt  least  squares  fit23  and  the  wavelength 
shift  of  the  two  centers  yields  the  axial  (velocity, 
while  the  width  of  each  Gaussian  yields  the 
temperature.  A  number  of  possible  uncertainties 
enter  into  this  analysis.  A  focused  pulsed  laser 
beam  creates  a  very  high  instantaneous  energy 
density,  which  is  needed  to  excite  the  two-photon 
transition,  but  also  may  induce  other  non-linear 
processes!  Multi-photon  ionization  (MPI)  will 
occur,  and  partial  saturation  of  the  transition 
becomes  possible,  even  though  the  two-photon 
absorption  cross  section  is  very  small.  In  addition, 
amplified  spontaneous  emission  (ASE)  will  take 
place,  since  a  population  inversion  is  created 
between  the  pumped  n=3  level  and  the  nearly  empty 
n-2  level.  The  population  inversion  can  cause  gain 
to  occur  for  any  photon  emitted  in  the  forward  or 
backward  direction  along  the  laser  beam.24'25,26  If 
ASE  or  MPI  become  large,  they  become  a  non- 
negligible  loss  mechanism  for  the  LIF  process. 
Saturation  will  cause  power  broadening  of  the 
Doppler  profile  of  the  transition,  and  saturation, 
ASE  or  MPI  will  cause  the  power  dependence  of  the 
LIF  signal  to  deviate  from  the  expected  behavior 
proportional  to  power  squared. 

We  have  performed  a  power-dependence  study, 
which  is  summarized  in  Figure  6  and  Figure  7.  We 
find  that  when  the  laser  beam  to  the  arcjet  is 
attenuated  to  below  0.03  mJ  per  pulse  for  the  radial 
beam,  the  LIF  is  proportional  to  the  power  squared, 
as  indicated  by  the  slope  of  2  on  a  log-log  plot.  At 
these  lower  laser  powers,  the  ASE  signal  has  fallen 
away  drastically  from  its  large  size  and  does  not 
appear  to  attenuate  the  LIF  signal.  At  higher  laser 
powers  reduction  in  LIF  signal  from  the  power 
squared  proportionality  is  observed.  In  addition, 
the  temperatures  measured  (shown  in  Figure  7) 
appear  to  increase  at  higher  laser  powers  where 
saturation  broadening  becomes  a  significant 


component  of  the  linewidth.  At  lower  powers,  the 
temperature  remains  constant  at  its  true  value. 

Figure  8  shows  the  differences  in  the  lineshape 
between  the  measured  ASE  signal  and  the  measured 
LIF  signal  in  wavelength  space.  Due  to  these 
differences,  compensating  for  the  ASE  loss 
mechanism  through  measurement  of  the  ASE  was 
deemed  unreliable.  Efforts  were  made  to  operate  at 
laser  powers  where  ASE  loss  was  insignificant  and 
did  not  cause  deviation  from  the  power  squared 
behavior  in  the  LIF  signal.  The  temperature  and 
number  density  data  presented  here  are  all  obtained 
with  a  laser  energy  of  about  0.03  mJ  per  pulse. 
Laser  powers  lower  than  this  were  found  to 
approach  signal  to  noise  limits  for  our  configuration, 
most  notably  when  probing  away  from  the  center  of 
the  arcjet  nozzle. 

Another  broadening  mechanism  that  is  present  is 
Stark  broadening  of  the  transition  due  to  the  free 
electrons  present  in  the  arcjet  plume.  We  do  not 
have  any  direct  measurement  of  the  profile  of  the 
electron  number  density  for  our  arcjet,  but  it  has 
been  measured  for  nearly  identical  conditions22  to  be 
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less  than  2  x  10  cm  at  maximum  in  the  center  of 
the  nozzle  exit  plane.  This  value  of  ne  would  cause 
a  Stark  width  of  0.002A  for  the  L|3  line.27  A 
simulated  Voigt  profile  using  this  Lorentzian  width 
shows  that,  for  a  typical  measured  linewidth  of 
0.029A,  accounting  for  the  Stark  broadening  would 
cause  the  temperature  from  the  Doppler  portion  of 
the  linewidth  to  go  down  from  1600K  to  1490K. 
This  represents  the  likely  maximum  uncertainty  in 
our  temperature  due  to  Stark  broadening,  and  will  be 
greatest  at  the  center  of  the  nozzle  where  ne  is 
largest. 

A  small  amount  of  line  broadening  is  caused  by  the 
linewidth  of  the  dye  laser  beam  which,  unlike  that  of 
cw  ring  dye  lasers,  is  non-negligible  compared  with 
the  atomic  linewidth.  The  laser  linewidth  at  205  nm 
has  been  measured  by  taking  an  LIF  spectrum  of  NO 
gas  in  the  cell.  Since  NO  is  much  heavier  than 
hydrogen,  and  is  at  room  temperature,  the  width  of 
the  rotational  lines  in  the  LIF  spectrum  reflect  only 
the  linewidth  of  the  laser  (as  long  as  care  is  taken  to 
use  unblended  rotational  lines).  This  has  shown  that 
the  laser  has  a  width  of  0.28  cm’1  or  0.012  A  at  205 
nm.  The  effect  of  the  laser  width  is  taken  into 
account  using  the  following  relation:28 
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(Av  )  =  (Avd  +  2Av,)2 

where  AvD  is  the  true  Doppler  width  of  the  line  (in 
cm'1),  Av,  is  the  width  of  the  laser,  Av  is  the 
measured  transition  width,  and  the  factor  of  two  is 
due  to  the  two-photon  probe  method. 

Velocity  is  measured  by  observing  the  Doppler  shift 
of  the  absorption  line  and  determining  the  velocity 
relative  to  the  incoming  photons  responsible  for  the 


where  vis  the  velocity,  AX  is  the  wavelength  shift, 
Xq  is  the  zero  velocity  line  center,  and  c  is  the  speed 
of  light.  Determining  the  line  center  of  a  profile 
with  the  fitting  process  discussed  earlier  is  relatively 
simple  and  quite  accurate.  The  uncertainty  from  the 
data  analysis  is  less  than  the  symbol  size  of  the  data 
shown  in  the  following  velocity  figures. 

Figure  9  shows  the  radial  distribution  of  axial 
velocities  0.4mm  from  the  arcjet  nozzle  exit.  The 
velocities  were  measured  with  a  laser  beam  entering 
the  plume  of  the  arcjet  axially.  Figure  10  shows  the 
corresponding  distribution  of  the  radial  component 
of  velocity  when  the  beam  is  brought  to  the  plume 
radially.  Data  for  the  plot  were  taken  over  several 
days  and  the’  scatter  may  indicate  day  to  day 
operational  changes  in  the  arcjet  as  well  as 
measurement  and  alignment  uncertainties.  It  should 
be  noted  that  velocity  profiles  of  the  ground  state 
hydrogen  atoms  are  in  close  agreement  with 
previous  measurements  using  excited  state  LIF.i5 

The  profiles  are  extremely  sensitive  to  centering,  in 
both  vertical  and  horizontal  directions  and  an 
alignment  procedure  is  performed  prior  to 
evacuation  of  the  chamber.  After  arcjet  operation 
commences  and  a  period  of  heating  (approximately 
30  minutes  to  an  hour)  takes  place,  sample  velocity 
distributions  are  taken  horizontally  and  vertically  to 
center  the  optical  detection  with  the  peak  velocity  of 
the  arcjet  plume.  This  usually  requires  an 
adjustment  of  0.25  mm  to  0.5  mm  from  the  initial 
alignment  at  atmosphere  in  one  or  both  of  the  axes 
that  make  up  the  nozzle  exit  plane.  The  y  (axial) 
axis  zero  is  defined  to  be  where  the  arcjet  nozzle 
moves  into  the  laser  beam  and  reduces  the  LIF 
signal  in  half.  Since  the  laser  is  focused  to  a  spot 


size  measured  to  be  .040  mm  ±0.01  mm,  this 
position  can  be  set  with  an  uncertainty  of  at  least 
±0.1  mm.  This  is  checked  periodically  to  insure  that 
it  remains  constant.  During  the  initial  wanning  up 
period  of  the  arcjet,  approximately  0.5  mm  of 
expansion  is  observed  in  the  1  kW  arcjet  nozzle. 

Figure  11  shows  a  temperature  profile  at  0.4  mm 
from  nozzle  exit.  The  results  indicate  a  peak 
temperature  around  1600K  -  1800K  in  the  center 
dropping  to  1000K  -  1200K  over  the  5  mm  radius  of 
the  nozzle  exit.  The  low  signal  to  noise  of  the 
fluorescence  (when  low  laser  power  is  used,  as 
described  above)  makes  measurement  of  the 
Doppler  width  difficult,  especially  at  the  edges, 
though  Doppler  shift  measurements  are  less 
affected.  Temperature  is  based  upon  the  square  of 
the  measured  Doppler  width  and  can  be  written  in 
terms  of  width  in  the  following  manner:29 

5X,  =  7.l6xlO“7/l 

where  AX  is  the  linewidth  due  to  Doppler 
broadening,  Xo  is  the  line  center,  M  is  the  molar 
mass  of  the  gas,  and  T  is  the  translational 
temperature.  The  constant  in  the  equation  is  based 
upon  the  values  of  the  speed  of  light  and  the 
molecular  gas  constant. 

Substantially  greater  uncertainty  is  associated  with 
temperature  than  with  velocity.  The  uncertainty  due 
to  scatter  can  be  seen  in  the  plot,  which  includes 
data  from  several  different  days. 

Figure  12  shows  our  observed  atomic  ground  state 
temperature  measurements  with  values  taken  from 
published  excited  state  temperature  measurements30 
and  molecular  rotational  temperatures."  While  the 
temperatures  predicted  by  the  ground  state  hydrogen 
atoms  are  quite  similar  to  the  rotational  molecular 
temperatures,  our  temperature  measurements  of  the 
ground  state  hydrogen  atoms  are  about  a  factor  of 
two  cooler  in  the  center  than  those  of  the  excited- 
state  hydrogen  atoms,  possibly  indicating  a 
significant  non-equilibrium  between  the  excited 
state  and  ground' state  hydrogen  atoms. 

The  dye  laser  was  tuned  to  656  nm  in  an  attempt  to 
directly  measure  the  excited  state  hydrogen  atoms  in 
our  arcjet  and  ensure  that  the  apparent  non¬ 
equilibrium  is  real  and  not  due  to  differences 
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between  out  system  and  Stanford's.  The  Ha  LIF 
was  detected,  but  the  power  broadening  due  to  the 
focused,  pulsed  laser  was  very  severe.  Decreasing 
the  laser  power  enough  to  alleviate  the  power 
broadening  led  to  an  unacceptable  signal-to-noise 
ratio,  so  this  measurement  could  not  be  pursued. 

Figure  13  shows  two  sample  fluorescence  decay 
curves,  one  measured  at  the  center  of  the  nozzle  exit 
plane,  and  the  other  measured  near  the  edge  of  the 
nozzle  exit  plane.  Change  in  the  lifetime  represents 
the  effect  of  collisional  quenching  in  the  3-2 
transition.  The  curves  are  obtained  by  digitizing  the 
LIF  signal  from  the  PMT  using  a  500  MHz 
oscilloscope.  Also  shown  in  the  figure  is  a  curve 
from  the  laser  scatter  only  (with  optical  filters 
removed  from  the  PMT)  to  show  the  time  resolution 
of  the  system  with  no  LIF  present.  Decay  curves 
such  as  Figure  13  are  analyzed  by  starting  at  10  ns 
after  the  laser  pulse  and  fitting  the  signal  to  an 
exponential  decay. 

Figure  14  shows  results  for  lifetime  measurements 
from  several  separate  days  as  a  function  of  position 
along  the  nozzle  exit.  The  lifetimes  measured  in  the 
arcjet  plume  indicate  the  collisional  quenching  due 
to  electrons,  H2,  or  most  likely  some  combination  of 
both.  The  LIF  signal  is  only  proportional  to  number 
density  if  it  is  corrected  for  losses  due  to  quenching. 

The  need  for  correcting  for  collisional  quenching  is 
evident  in  Figure  15  where  the  density  prediction 
before  and  after 'the  correction  for  quenching  is 
seen.  In  this  case  the  correction  is  significant  in  the 
center  region,  and  minimal  at  the  edge  of  the  profile. 

Figure  16  shows  the  calibrated  number  density 
profile  of  hydrogen  atoms  near  the  exit  plane  of  the 
nozzle  (0.4mm  from  nozzle  exit  plane).  This  is 
obtained  from  the  area  under  the  curve  of  the 
hydrogen  atom  profile  (at  low  laser  power  to  avoid 
ASE,  MPI  and  saturation),  and  corrected  by  the 
fluorescence  yield  at  each  location.  Separately 
calibrated  data  sets  from  many  different  days  are 
included  in  the  plot.  The  peak  of  the  density  profile 
is  at  the  center  of  the  nozzle  exit  plane  and  indicates 
a  number  density  of  1  x  1016  cm'3  and  falls  off  about 
one  and  a  half  orders  of  magnitude  over  the  5  mm  to 
the  nozzle  exit  wall. 

Also  shown  in  the  figure  are  the  predicted  exit  plane 
number  density  distributions  from  three  arcjet 


models.  The  Boyd  model  is  authored  by  I.  D.  Boyd 
at  Cornell  University.9  It  is  a  Direct  Simulation 
Monte  Carlo  (DSMC)  particle  code  that  predicts 
physical  properties  through  most  of  the  arcjet  nozzle 
and  out  into  the  plume.  The  Butler  model  is 
authored  by  G.  ,\V.  Butler  at  01  in  Aerospace 
Corporation  and  is  a  Navier  Stokes  (NS)  code  that 
computes  the  flow  properties  throughout  the  arcjet 
nozzle,  stopping  at  the  nozzle  exit.5''1'  A  separate 
NS  code  written  by  Megli,  Krier,  and  Burton  from 
the  University  of  Illinois  predicts  results  that  are 
also  shown  in  the  figure.'12  Like  the  Butler  model 
the  calculations  extend  to  the  nozzle  exit,  but  not 
into  the  plume. 

The  number  density  differences  between  the  data 
and  the  model  indicate  almost  an  order  of  magnitude 
higher  atomic  density  than  predicted  by  the  two 
Navier  Stokes  models  and  a  factor  of  three  higher 
than  the  DSMC  prediction. 

Shown  in  Table  1  is  a  comparison  of  the  predictions 
for  both  atomic  hydrogen  and  molecular  hydrogen 
densities  for  each  of  the  models.  In  addition 
experimental  results  on  species  density  from  this 
work  and  that  of  Beattie  and  Cappelli  at  Stanford 
University  are  shown.11 

Note  that  the  predictions  for  all  three  models  agree 
well  with  each  other  and  with  experiment  for 
molecular  hydrogen  density  at  the  centerline  of 
nozzle  exit.  As  seen  in  the  previous  figure,  the 
differences  in  atomic  number  density  can  be  quite 
large.  These  atomic  discrepancies,  in  light  of  the 
molecular  density  agreements,  indicate  different 
molecular  dissociation  fractions  predicted  by  the 
codes  and  by  the  combination  of  the  two 
experiments  at  the  center  of  the  nozzle  exit. 

The  two  experimental  efforts  indicate  a  dissociation 
fraction  as  large  as  50%  at  the  nozzle  exit  center. 
This  is  an  approximation  for  dissociation  as  the 
error  bars  on  both  experiments  are  quite  large  and 
the  experiments  were  done  on  two  different  arcjets, 
though  every  effort  was  made  to  operate  at  similar 
arcjet  conditions. 

The  predicted  dissociation  fraction  from  each  of  the 
models  is  substantially  less  ranging  from  5%  to  18% 
depending  upon  the  model.  Accurately  predicting 
the  dissociation  fraction  could  be  closely  related  to 
understanding  the  amount  of  frozen  flow  energy  loss 
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in  the  arcjet  thruster  and  an  important  step  in 
improving  arcjet  efficiency  from  the  current  30%  to 
35%  range. 

Figure  16  shows  a  two  dimensional  profile  of  the 
hydrogen  atomic  number  density,  data  taken  at  .4 
mm  downstream  of  nozzle  exit.  This  figure  shows  a 
generally  axisymmetric  flow  with  non-smooth 
features  at  the  edges  of  each  density  range.  The 
rough  contours  are  caused  partially  by  the  1/distance 
algorithm  used  to  create  the  contours  from  a 
Cartesian  data  set  that  was  large  and  intensive  to 
take,  but  still  less  in  quantity  of  points  than  perhaps 
desirable  for  identifying  precise  gradients.  In 
addition,  the  noise  in  the  data  as  seen  in  the  earlier 
density  plot  of  Figure  16  is  significant  and  can  cause 
the  appearance  of  shifting  gradients  in  this  type  of 
two  dimensional  representation. 

Figure  17  shows  the  same  data  in  the  two 
dimensional  plot  all  collapsed  onto  a  single  axis 
with  radius  as  the  ordinate.  Each  line  is  a  single 
"spoke"  of  the  two  dimensional  data  shov/n  in 
Figure  17.  The  lines  are  shown  together  to  illustrate 
that  the  spread  between  the  different  radial  profiles 
is  within  the  scatter  that  is  seen  when  examining  just 
one  radial  profile  on  several  different  days.  This 
indicates  that  the  flow  appears  to  be  quite 
axisymmetric  within  the  experimental  uncertainty. 

Examining  the  behavior  of  the  density  and 
temperature  as  the  flow  exits  the  nozzle,  Figure  19 
shows  the  density  and  temperature  gradients  from 
0.2mm  past  nozzle  exit  to  30  mm  downstream  along 
the  nozzle  centerline.  One  of  the  benefits  of  looking 
at  these  trends  was  to  examine  the  validity  of 
comparing  data  taken  at  0.4  mm  downstream  to  the 
properties  right  at  nozzle  exit.  Since  the  optical  data 
for  this  experiment  could  not  be  observed  inside  the 
nozzle,  many  data  points  were  taken  close  to  the 
nozzle  exit  and  the  axial  position  axis  is  shown  on  a 
logarithmic  scale.  Note  that  while  the  density  does 
appear  to  decrease  as  it  exits,  the  difference  between 
the  data  taken  at  0.4  mm  and  what  would  be 
extrapolated  back  to  nozzle  exit  is  well  within  the 
error  bars  shown  in  Figure  16,  and  one  expects  that 
the  data  would  not  appear  significantly  different  if 
the  exact  nozzle  exit  was  able  to  be  probed  with  this 
technique.  In  addition,  if  the  data  fails  to  predict  the 
correct  atomic  number  density  right  at  the  nozzle 
exit,  it  would  appear  from  the  trend  in  this  figure  to 
underpredict  the  number  of  hydrogen  atoms  at 


nozzle  exit  leading  to  an  even  greater  discrepancy 
with  the  predictions  from  the  computational  models. 

Conclusions 

Using  two-photon  LIF,  hydrogen  ground  state 
number  densities,  velocities  and  temperatures  have 
been  measured  near  the  arcjet  nozzle  exit.  General 
agreement  with  previous  nozzle  exit  velocity 
distributions  is  seen,  with  no  effect  of  translational 
slip  being  observed.  Temperature  measurements  of 
the  atomic  ground  state  appear  to  be  significantly 
cooler  than  those  previously  measured  in  excited 
states  indicating  non-equilibrium  behavior. 

Number  density  measurements  at  the  nozzle  exit  are 
made,  carefully  avoiding  such  loss  mechanisms  as 
multi-photon  ionization,  amplified  spontaneous 
emission,  and  laser  saturation.  Fluorescence 
quenching  was  measured  and  corrected  for. 

The  number  density  distributions  are  compared  with 
computational  modeling  data  and  indicate  a 
significantly  higher,  more  peaked  number  density 
than  predicted.  Including  recent  molecular  density 
data,  a  higher  dissociation  fraction  is  observed  than 
is  currently  predicted  by  the  computational  models 
indicating  that  frozen  flow  energy  loss  may  be  a 
more  significant  contribution  to  the  arcjet's  low 
efficiency. 

Axial  symmetry  of  the  density  at  nozzle  exit  is 
observed  within  the  scatter  of  the  measurement  and 
both  density  and  temperature  axial  trends  are 
observed  indicating  that  data  taken  at  0.4mm  is 
likely  representative  of  properties  right  at  the  nozzle 
exit  and  data  across  the  nozzle  exit  in  one  direction 
is  representative  of  the  entire  plane. 
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Figure  i.  Oc termination  of  atomic  hydrogen  density  and  emitted  signal  throtieh  titration  of  NOjinto  calibration 


FigureS.  Sample  2PLIF  spectrum  of  the  hydrogen  atom,  showing  the  calibration  signal  from  the  discharge  cell  and 
ih«  si(«na>  from  (he  arcjct  plume. 
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Figure  9.  Profile  across  oottie  exil  plane  (0.4  mm  from  exit)  of  the  axial  velocity  component  of  the  ground  *,3,c  ' 
hydrogeo  atoms  containing  data  from  several  diffcscnt  days.  For  comparison,  corresponding  data  (same  arcjet 
conditions,  same  profile  location)  for  LIF  of  the  electronically  excited  hydrogen  atoms  from  Stanford  University.** 


Figure  d.  Arcjct  study  of  UF  signal  dependence  on  laser  power  compared  to  ASE  power  dependence,  which  has  a 
much  thurptr  response  to  laser  power. 
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Figure  10.  Profile  across  no  ale  exit  plane  (0.4  mra  from  exit)  of  the  radial  velocity  component  of  the  ground  state 
hydrogen  atoms  containing  data  from  several  different  days  where  the  radial  position  is  measured  parallel  to  the 
incoming  laser  beam. 


Figure  7.  Apparent  temperature  for  the  arcjet  at  changing  laser  powers 


Laser  Wavelength  (ran) 

Figure  8.  Narrower  ASE  spectrum  compared  with  LIF  spectrum. 


Figure  II.  Profile  across  norzle  exit  plane  (0.4  mm  from  exit)  or  the  translational  temperature  or  the  ground  state 
hydrogen  atoms  containing  data  from  two  different  days.  The  error  bars  arc  longer  on  the  lower  ride  to  reflect  the 
estimated  amount  of  Stark  broadening. 
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Figure  12.  Temperature  data  shown  in  Figure  II,  with  corresponding  daia  (same  arcjct  conditions,  same  profile 
location)  for  LIF  of  the  electronically  excited  hydrogen  atoms''4  and  Raman  molecular  rotational  temperatures, 
both  from  Stanford  University, 


